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We investigate the role played by electron traps on adiabatic charge transport for graphene 
nanoribbons in the presence of an acoustically induced longitudinal surface acoustic wave (SAW) 
potential. Due to the weak longitudinal SAW-induced potential as well as the strong transverse 
confinement by a nanoribbon, minibandsof sliding tunnel-coupled quantum dots are formed so that 
by varying the chemical potential to pass through the minigaps, quantized adiabatic charge trans- 
port may be obtained. We analyze the way that the minigaps may be closed, thereby destroying 
the likelihood of current quantization in a nanoribbon. We present numerical calculations showing 
the effects due to electron traps which lead to localized-trap energy levels within the minigaps. 
Additionally, for comparison, we present results for the minibands of a corrugated nanoribbon in 
the absence of a SAW. 

I. INTRODUCTION 

A considerable amount of research work has been carried out so far on the design and improvement of electronic 
devices which are based on the use of quantized adiabatic charge transport. P^10| Moreover, under a surface-acoustic 
wave (SAW), the inelastic capture and tunneling escape effects on the non-adiabatic transport of photo-excited charges 
in quantum wells was also investigated. |11) The underlying challenge is to produce a device with an accuracy for the 
quantized current of one part in 10* on the plateaus. When this goal is achieved, one application of this device would 
be in metrology for standardizing the unit of current. 

At the present time, a SAW is launched on a piezoelectric heterostructure, such as GaAs/AlGaAs, and GHz 
single/few-electron pumps have been gaining close scrutiny due to the fact that the measured currents lie within the 
nanoamp range, high enough for the measured current to be suitable as a current standard. However, these pumps 
have so far been capable of delivering electrons/holes in each cycle of a sliding dynamic quantum dot (QD), giving 
rise to a quantized current with an accuracy of one part in 10^ as reported in Refs. [3HH]. Interestingly, in Ref. ^1^ 
a measurement was carried out of the noise accompanying a 3-GHz SAW pump. It was observed in this experiment 
that the current near the lowest plateau, corresponding to the transfer of one electron per SAW cycle, is dominated 
by shot noise. However, away from the plateau, the noise is attributed to electron traps in the material. There have 
been some attempts to increase the flatness of the plateaus by applying magnetic fields. [T31 E] 

Some time ago, a proposal was put forward by Thouless [10] which would make use of quantized adiabatic charge 
transport. This adiabatic approach involves the use of a one-dimensional (ID) electron system subjected to a slowly- 
sliding periodic potential. Relatively simple analysis indicates that in such a ID system minigaps are generated in 
instantaneous electronic spectra as a function of the SAW amplitude. With the use of a gate, the chemical potential 
can be varied by applying a voltage to the gate. Consequently, when the chemical potential lies within a minigap, there 
will be an integral multiple of electron charge transported across the system during a single time period. [T] In other 
words, by combining with the strong transverse confinement of a nanoribbon, the weal longitudinal SAW potential has 
induced a series of dynamic (sliding) tunnel-coupled QDs whose impenetrable "wall is constructed through destructive 
interference of the electronic wave functions around a minimum of the SAW potential. In principle, such an adiabatic- 
transport device could provide an important application, like a current standard. Talyanskii, et al. [T] investigated the 
physical mechanisms of quantized adiabatic charge transport in carbon nanotubes for a SAW to produce a periodic 
potential required for miniband/minigap formation. 

In the presence of a SAW, the scattering effects from impurities embedded in a ID electronic system are expected to 
play an important role on the flatness of a current plateau. The current quantization should be completely smeared out 
when the level broadening from impurity scattering becomes comparable to the minigaps of dynamic tunnel-coupled 
QDs. On the other hand, we can also simulate localized electron traps by superposing a series of negative ^-potentials 
onto a SAW potential within each spatial period. Consequently, we expect a set of localized trap states occurring 
within the minigaps of dynamic QDs. This provides an escape channel for the QD-confined electrons being carried 
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by the SAW. This trap mechanism is quite different from the impurity one [T* where a spatial average with respect to 
the distribution of impurities within a dynamic QD is inevitable due to a SAW. 

In this paper, we consider a ID Dirac-like electron gas in a graphene nanoribbon in the presence of a SAW. We will 
introduce two mechanisms for miniband formation. First, the nanoribbon is modulated by a longitudinal potential 
from a SAW. Secondly, the nanoribbon is periodically corrugated. We notice that the second mechanism does not 
lead to a quantized current but instead produces traps for Dirac electrons, thereby limiting electron mobility. Our 
numerical calculations reveal that localized electron trap states are an effective mechanism to adversely affect the 
adiabatic transport because the localized-trap levels lying within the minigaps are very sensitive to the phase of either 
the SAW or the corrugation-induced potential. Varying the weight or the position of the ^-potential leads to different 
positions of localized trap levels within the minigaps of the nanoribbon. Therefore, these inevitable fluctuations of 
the trap potential in a realistic system would most likely impede the current quantization. 

The rest of the paper is organized as follows. In Secjllj we present the formalism for calculating band structure 
with localized trap states for nanoribbons in the presence of a SAW. In Sec. |III[ numerical results for nanoribbons 
in the absence/presence of a SAW and those for corrugated nanoribbons in the absence of a SAW are presented to 
demonstrate and explain the localized trap states within the minigaps. The conclusions drawn from these results are 



briefly summarized in Sec. IV 



II. MINIBAND STRUCTURE WITH LOCALIZED TRAP STATES 



The work done by Talyanskii, et al. pL^ on quantum adiabatic charge transport focused on the coupling between a 
semimetallic carbon nanotube and a SAW. The electron backscattering from the SAW potential is used to induce a 
miniband spectrum. The electron interactions enhance the minigaps thereby improving current quantization. The 
effect due to impurities in the carbon nanotube is averaged by a SAW potential. 

For the cases of a semimetallic carbon nanotube, semiconducting carbon nanoribbon with applied SAW potential 
and corrugated nanoribbon, the energy levels are given by the spectra of discretized ID Dirac Hamiltonian (see the 
Appendix A for detailed derivations) 
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The eigenvalue problem is defined within the spatial interval < x < 2'K/k and assumes periodic boundary conditions. 
In this notation, k stands for either the wave number ksAW of the SAW potential or the wave number kc of an effective 
potential induced by the corrugation. The discretization of the Hamiltonian is provided by the mesh a;„ = nSx with 
71 e . . . iV - 1 and (Sx = 2Tr/kN. 

For either the carbon nanotube or graphene nanoribbon, the parameters for the Hamiltonian matrix in Eq. ([l]) are 
given by 



a„ - , (2) 
ihvp 



where we have introduced the SAW and impurities combined phase 
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a{x) = y [VsAw{u) ~ Vtrap{u)] = Xcos{ksAwx) + Vo9{x - xo) , (3) 



with A = 2A/{hvFksAw) the normaUzed SAW amphtude, vp is the Fermi velocity of graphene. Additionally, 
Vq — Vtrapl {fivpksAw) and x^) denote the normalized trap-potential amplitude and position of the short-range 
dynamic trap for electrons, respectively, and the trap is sliding together with the SAW potential. The mass term 
involving A is the original energy gap for the system in the absence of a SAW. In case of a nanotube, A may be 
generated by a magnetic field. For a nanoribbon, the gap is structural for the semiconducting nanoribbon hvpky"'^ 
with fcy™^ being the transverse electron wave number due to finite size across the ribbon. For nanoribbons, the explicit 
form for the phase introduced in Eq. ([3| can be found from the Appendix A. 

As far as the minigaps are concerned, the effect due to the SAW potential on the nanoribbon may be compared 
with corrugation. A sinusoidal corrugated semiconducting ribbon can be mapped on to a flat ribbon. The mapping 
introduces an additional cti term into the Dirac equation, as described in Appendix A. This yields 



hvj. 



-V^C^kl sin(2fcc x„ 



A,(a;)[2-FC2fc2 cos(2fc,x„)-hC2fc2]3/2 4A3(a;) 

ihvp 



C k^ • /n 7 \ 

I . . „ , , sm(2fccX„) 



(4) 



&„ = ihipk'f^ e-2»"<=(^") , 



Ac(a;„) 25x 



where the corrugation-induced gap is given by Ac(x„) = + C^k^ cos2(fcca:„), with C being the normalized amplitude 
of the corrugation, and fc^™'' is the quantized wave number across the nanoribbon. 



III. NUMERICAL SIMULATION AND DISCUSSIONS 



In our numerical calculations, all the energies in Figs. [T] and [2j such as e, A and are normalized to hvpksAW ■ 
The SAW potential amplitude A is also normalized to hvpksAW- Additionally, all the energies in Fig.[3j such as e and 
A, are measured in units of hvpkc, and the corrugation amplitude C is normalized to l/kc- Besides, the transverse 
wave number fcy™^ in all the plots is scaled by 27r/3ao. In this way, we are able to draw some universal conclusions 
concerning the effects due to minigaps. 

In Fig.[l} we compare the energy band structure of nanoribbons for two values of A in the absence of electron traps. 
Two values of ky"^^ were chosen (light and dark colors) to describe the two lowest energy levels (see the discussions 
in the Appendix A). The minigaps are generated by a sliding dynamic QD and they oscillate as a function of the 
SAW amplitude as may be verified using perturbation theory, vanishing at values close, but generally not equal, 
to the roots of Bessel functions. Increasing or decreasing the value of A results in a shift of the nodes on the graph 
as evidenced by comparing our results in Fig.jl] Therefore, A determines not only the magnitude of the original gap 
in the absence of a SAW but also the size of the minigaps in the presence of a SAW. Higher energy minigaps are 
partially closed by the energy levels corresponding to a larger value of /cy™"* (not shown here). 

We now introduce electron traps into our nanoribbon by superposing a negative (5-potential onto the SAW potential 
so as to simulate a short-range Coulomb interaction. In this case, the position of the trap is fixed in the moving 
SAW frame of reference, which is quite different from embedded impurities in a nanostructure. In the moving SAW 
frame, the embedded impurities are moving against the dynamic QDs created by both the transverse dimension of 
the nanostructure and the longitudinal SAW potential. This results in an average of the impurity effects with respect 
to these dynamic QDs in the longitudinal direction. As seen in the results presented in Fig. [2] localized trap states 
occur within the minigaps once the weight of the trap potential Vq becomes larger than A/{hvpksj^yir). Relative 
energy value of these trap states in the presence of SAW is sensitive to the position xa of the trap within a dynamic 
QD. If we set A = Vq, then the contribution to a{x) from the trap located in the nodes of the SAW potential is fully 
compensated by the cosine term in Eq. ([3]). As a result of this compensation, the localized trap states will disappear 
from the gap and minigap regions. If we extend the single-trap model employed in this paper to a uniform distribution 
of traps, the fiuctuations in the phase term a{x) [see Eq. p|] would fill up the entire minigap region with a delocalized 
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trap band. Consequently, the adiabatic approximation may not be applicable. In other words, to satisfy the adiabatic 
assumption, one must have dominance of the SAW potential, i.e., A 3> Vq must be satisfied. 

We compare the results for S AW-based dynamic QDs in Figs.[T] and [2] with those for static QDs created by corru- 
gation on a graphene nanoribbon in the absence of a SAW and electron traps. This we do by displaying in Fig. [3] the 
minigaps induced by the corrugation. We find from the figure that minigaps only exist for finite but small values of 
the corrugation amplitude C. This means that these minigaps are generally much less than those induced by a SAW. 
As a matter of fact, the existence of non- vanishing diagonal terms a„ given in Eq. Q, effectively mitigates the phase 
fluctuations in the off-diagonal terms &„. This keeps the minigaps open and the energy spectra robust even after traps 
have been introduced to cause a fluctuation in the phase term ac{x). 

Finally, let us assume that a narrow channel is formed within a two-dimensional electron-gas layer lying in the xy- 
plane. We will neglect the finite thickness of the quantum well for the heterostructure in the z-direction and consider 
the electron motion as strictly two dimensional. We will employ one of the simplest models for the gate-induced or 
etched [15] confining electrostatic potential. In this way, a ID channel is formed on the two-dimensional electron-gas 
layer and the dynamics of massive electrons can be modeled by a discretized ID Schrodinger equation. However, from 
numerical results (not shown here), we find no evidence of the minigaps for this model, i.e., the minigaps are the 
characteristics of Dirac fermions. 



IV. CONCLUDING REMARKS 

In conclusion, we have calculated in this paper the energy band structure for graphene nanoribbons, embedded 
with a single electron trap, upon which a SAW is launched. Our results show that localized trap states appear in 
the minigaps. More importantly, the location of the trap state-energy level is determined by the positions of the trap 
with respect to the phase of the sinusoidal SAW. Consequently, the adiabatic approximation might not be appropriate 
whenever the minigap is less or comparable with the weight of a short-range 5-potential for the trap (see Fig. [2]). On 
the other hand, in Fig.[T] where there are no electron traps, a larger value of A in the energy spectrum leads to a 
substantial increase in the number of minigaps as well as the ballistic current quantization. Periodic corrugation of 
the nanoribbon may be used instead of a SAW as a mechanism for inducing minigaps. Those are expected to be less 
sensitive to the presence of charged impurities or electron trap potentials. 
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Appendix A: Energy Band Calculations in Absence of Electron Traps 

1. Carbon Nanotubes 

The electron eigenstates in a semi-metallic nanotube are described by a ID Dirac equation. For simplicity, a 
noninteracting system is considered here. Under the stationary approximation, the single particle energy spectrum 
e(fc) is obtained from the following perturbed ID Dirac equation 



e(k) tpaix) = —ihvp — ^ h A-!/)fl(a;) + Asm{kx) ipaM , 

ox 

e{k)^fi{x) = ihvp ^'^^^^^ + A-0Q(a;) + Asin(fca;) i/'/sla:) . 



(Al) 
(A2) 



In this notation, k represents the electron wave number along the nanotube, vp is the Fermi velocity of Dirac electrons, 
A is the SAW amplitude, A is the energy gap of the system in the absence of a SAW, a and /? label the two sublattices 
of graphene from which the nanotube is rolled. In addition, we require k = ksAW to satisfy momentum conservation, 
where ksAW is the wave number of a SAW propagating along the nanotube. To explore the miniband structure due 
to quantum confinement in the radial direction, a gauge transformation is implemented and is defined by 



ll}fi(x) 



g(»/2)Acos(fca;) ^^(-^) 



(A3) 



where A — 2A/{hvFk) is the normalized SAW amplitude. Substituting Eq. (A3) into Eqs. (Al) and (A2), we obtain 



d 



+A e(-^/^)^^°<kx) _^ ^ sin(fca;) e^^/^)^ i/'^ {x) , 

:(fc) e(-V2)Acos(M ^ ^^^F^ 

+Ae('/2)^™^('=^) i,'p{x) + Asm{kx) e(-V2)Acos(M ^/^^^^ 



Ji/2)Acos(fcx) 



(A4) 
(A5) 



By introducing the following identities 



-ihv}. 



d_ 
dx 



Ji/2)\cosikx) 1 1 



-Asin(kx) Mx) - ihvF e^'/^)^™^^^^) ^^^M 

ox 

d ^M/2)Acos(My,^(a;)" ^ 



dx 



A ■ n \ ;' /- ^ I -fc- (-i/2)\cos(kx) (^"^'ai^) 

—Aiiin[kx)ip^[x) + invFe — ^ 



Eqs. (A5) and (Al) may be simplified as 

d^' (x) 

e{k) i,'{x) = -ihip -^-^ + A e-''^^-('=-) ^'^{x) 
^ ox 



e{k) ^P'^ix) = ihip "^^^ + A e'^^^C^-) ^'^{x) 



(A6) 



(A7) 



(A8) 
(A9) 



Furthermore, by employing the basis set ^'(x) = |'0^(a;), the above equations can be rewritten into a 

compact matrix form, given by 
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e{k)li''{x) ^ni/'ix) , (AlO) 

n = thvpffa ^+Aai e-'^^' ^^^^^^ . (All) 
ox 

We will solve the eigenvalue problem within the spatial interval < x < 2TT/k and introduce the iV-point mesh 
Xn = n Sx, where n Cz . . . N — 1 and Sx = 2TT/kN. In this way, the derivative on the mesh can be approximated by 



d'^' [xn) / dx — > [^'{xn+i) ~ {xn-i)]/ {'i.Sx) . Especially, on this spatial mesh, the Hamiltonian in Eq. (All) may be 
projected into the matrix given in Eq. ([T]). 

2. Graphene Nanoribbons 

Here, we consider an armchair graphene nanoribbon lying along the x-direction. The total number of carbon atoms 
(in both sublattices) across the ribbon is assumed to be M. The armchair edges mix up the graphene K and K' 
valleys so that the wave function becomes 

*™(a;, y) = e^'=«'"'^ *k^„(x) + 6"''=^'^ *k,^„(x) , (A12) 
where the transverse wave number is given by 

kl-^ = ^ + ^ , (A13) 

m ~ 0, ±1, ±2, • • • is an integer, L is the nanoribbon width, and ao = \/ia/2 {a w 1.42A) is the size of the unit cell 
in graphene. Nanoribbons with width L/ao = 3Af + 1 give rise to the following relation 

^ 2nm 2^ ^ ^ ( 2M + l + m \ 

y 6Mao + 3ao 3ao 3ao \ 2M + 1 J ' ^ ' 

and it is clear that the minimum energy occurs at ky = 0. Therefore, such nanoriibons are metallic. The 

next miniband corresponds to ky = (27r/3ao) [l/(2i\f + 1)]. On the other hand, for nanoriibons having width 
L/aQ — 3M (upper Eq.) or L/qq = 3M — 1 (lower Eq.), the two minimal values of fc^™'' are found to be 



4-2*'^) = (27r/3ao) [1/(6M + 1)] and 4"^*^"^^ = (27r/3ao) [-2/(6M + 1)] , 
4"'^^^ = (27r/3ao) [-l/(6Af - 1)] and 4"'^^+^' = (2^/3ao) [2/(6Af - 1)] . 

Those nanoribbons are semiconducting with the energy gap determined by 



(A15) 



3ao \6M ± 1 



The x-component of the wave function in Eq. (A12| may be determined by 



e{k)mKiK'),m{x) = n'^K{K')..m(x) . (^17) 

n = ±hvp f-iai ^ + ik^^^ aA+A sin(fca;) I , (A18) 



where the ± signs correspond to K and K' valleys. Additionally, the Hamiltonian in Eq. (Al8l can be transformed 



into the form in Eq. (All ) after applying the following unitary transformation 

U = (A19) 



^-ia(x) Qia{x) 
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where a{x) = —Acos(ksAwx)/{f^VFksAw)- As a result, the transformed Hamiltonian takes the form 



-id/dx -i4"^e*2a(x) 
i4"^e-*2"(=^) id/dx 



(A20) 



Formally, this Hamiltonian is equivalent to that in Eq. (All) after we applying the following substitutions 



2a{x) — > —Xcos{ksAwx) 



(A21) 



The valley sign ± does not change anything due to the mirror symmetry in the energy dispersion relation e{k) with 
respect to ±A:. 

3. Corrugated Nanoribbons 

We now turn to the case of a corrugated graphene nanoribbon whose modulation is sinusoidal with amplitude C 
and wavelength iir/kc along the x-axis. The model Hamiltonian for such a corrugated graphene nanoribbon has been 
given in Ref. [TB] as 



a) 7,- ' ^ i. ^ " df{x)/dx 



dx 



2 A,(a;) ' 



where 



A.(.0^/rTW(^ and i^(^) - ^^^^^^^ 
f{x) = C shi{kcx) . 
After applying the following unitary transformation 



U = 



72 [ 



and U'^ = 



72 



where 



(A22) 



ac{x) = — r-^ / du 



sin (2kcu) 



4 J [l+C2fc2cos2(fccM; 



i3/2 



_ 1 

~ 2 

the transformed Hamiltonian becomes 



V2 



72+C2fc2 cos(2fc,a:)+C2fc2 ^1+C2fc2 



ii' = hip 

where the complex function is defined by 



ik e~2i"c n*f^^ j_ '^/^^ 



^*( \ I '■a/ax 



(A23) 



G{x) 



'^/2C'^kl sin(2fcc X) 



Ae(a;) [2 + C2fc2 cos(2fc,x) + C2fc2]3/2 4A3(a:) 
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By making use of the finite-difference method for calculating d/dx along with the following basis set 



^ — {An-l, Bn-l,An, Bn, An+1, B„ + l} 

the transformed Hamiltonian matrix in Eq. ( |A23 1 may be projected as 







I 

A^{x„) 2Sx 







G* {Xn 


-l) 









2Sx 













n ? 

" A<:(x„_i) 2(5a; 

ifcj,e-2*«=(^") G*{Xn) 









Ae(2:„) 2Sx 



G(x„+i) 



(A24) 








i 

Ac(a;„) 25a; 
_jj[j.^g2iac(a;7i + l) 



The above Hamiltonian matrix has the same form as Eq. ([T]). 



A.(.^)25. *fc.e-2-^(-+^) G*K+i) 
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FIG. 1: (Color online) Scaled electron energy spectrum e{k) /(hvpksAw) in the absence of electron traps as a function of the 
normalized SAW amplitude A/ {fiVFksAw) for semiconducting nanoribbons subjected to an acoustically induced SAW potential. 
In this figure, we choose A/ {hvpksAw) = 1.0 (upper panel) and 1.5 (lower panel). Only the eigen-spectra arising from the 
two lowest dispersion curves in the absence of a SAW are displayed. Higher subbands contribute significantly at larger SAW 
amplitude. The lighter shaded regions arise from the lowest energy dispersion curves, whereas the darker shaded regions are 
associated with the energy dispersions of the next subband. 
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4/;iksAwi'F=l 0, £g/»sAwvjr=1.0 A/ftksAwVF=1.0, Eg/'iksAwi'f =1.0 




FIG. 2: (Color online) Scaled electron energy spectrum e{k) / [hvpksAw) in the presence of electron traps as a function 
of the normalized SAW amplitude A/ [hvpksAw) for semiconducting nanoribbons subjected to an acoustically induced SAW 
potential. In this figure, we chose the energy gap /S./ {hvpksAw) = 1-0 and the trap weight Vo = Vtrap/ {hvFksAw) ~ 1.0 for 
all four panels. As in Fig.[l] only the eigen-spectra due to the two lowest dispersion curves in the absence of the SAW are 
shown. Here, the position xoksAw / (Stt) of electron traps, moving with the SAW potential, is located at: 1/3 (upper-left panel), 
1/5 (lower-left panel), 1/7 (upper-right panel), and 1/11 (lower-right panel). Two induced trap-state energy levels within the 
energy gap are indicated by arrows in the lower-left panel for emphasis. 
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FIG. 3: (Color online) Minigap spectrum e{k) / (hvpkc), induced by a corrugation potential, as a function of the normalized 
modulation amplitude Ckc- In this figure, the two plots correspond to the two lowest values of the quantized transverse wave 
number ky"^^ . The graph at the top comes from the lowest quantized energy subband, whereas the graph at the bottom is due 
to the first-excited subband. 



